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When a ball with a concentric hole through it is suspended within a spin- 
ning tubular projectile the ball will rotate in such a way that its hole is 
aligned with the axis of spin of the projectile. The ball thus serves as an 
automatic obturator--plugging the projectile while the two bodies are within 
the gun barrel and rotating to open the tubular passage following exit from 
the muzzle. The motion of the ball is gyrodynamic in nature and highly depen- 
dent upon the external moments on the ball that arise because of its motion 
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relative to the spinning projectile. 

This report presents the results of an analytical and experimental 
study to provide an understanding of such motions in as complete a sense as 
is possible. The general equations of motion are formulated together with 
appropriate mathematical expressions for the external moments. These moments 
have their origins in the forces acting on the ball (gravity and aerodynamic 
forces are considered) which, in the presence of relative motion between 
the ball and the projectile, lead to sliding friction torques. Torques 
due to fluid shear are also considered but these are found to be relative- 
ly insignificant. 

Exact solutions are obtained using standard numerical techniques. In 
addition, a linear form is developed' and these solutions lead to useful 
approximations that are valid over a broad range of operating conditions. 

Experiments are described in which the general validity of the theoreti- 
cal models (exact and linear) is demonstrated for cases in which there is 
little or no aerodynamic load. The experiments, together with the linear 
approximation, lead to a semi-empirical method for determining the effective 
coefficient of sliding friction for such systems. 

A model is proposed to account for aerodynamic loading and sensitivity 
studies are conducted to determint the nature and scope of the influence 
of various design parameters upon systme performance. The linear approxi- 
mation is proposed as a useful design guide when applied with due awareness 
of its limitations. In addition, a design criterion is presented by means 
of which it is possible to avoid designs that lead to operation in a 
"hovering" region. Such operation leads to prolonged delays in ball opening 
time. Application of the design criterion leads, to ball/projectile 
designs that are optimum from the point of view of minimum ball opening time. 
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SUMMARY 



When a ball with a concentric hole through it is suspended within 
a spinning spherical cavity the ball will rotate in such a way that its 
hole is aligned with the axis of spin of the cavity. If it is mounted 
within a spinning tubular projectile, the ball can serve as an automatic 
obturator--plugging the projectile while the two bodies are within the 
gun barrel and rotating to open the tubular passage following exit from 
the muzzle. The motion of the ball is gyrodynamic in nature and highly 
dependent upon the external moments on the ball that arise because of its 
motion relative to the spinning cavity. 

This report presents the results of an analytical and experimental 
study to provide an understanding of such motions in as complete a sense 
as is possible. The general equations of motion are formulated together 
with appropriate mathematical expressions for the external moments. These 
moments have their origins in the forces acting on the ball (gravity and 
aerodynamic forces are considered) which, in the presence of relative 
motion between the ball and the projectile, lead to sliding friction torques. 
Torques due to fluid shear are also considered but these are found to be 
relatively insignificant. 

Exact solutions are obtained using standard numerical techniques. 

In addition, a linear form is developed and these solutions lead to use- 
ful approximations that are valid over a broad range of operating condi- 
tions . 

Experiments are described in which the general validity of the 
theoretical models (exact and approximate) is demonstrated for cases in 
which there is little or no aerodynamic load. The experiments, together 
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with the linear approximation, lead to a semi-empirical method for deter- 
mining the effective coefficient of sliding friction for such systems. 

A model is proposed to account for aerodynamic loading and sensi- 
tivity studies are conducted to determine the nature and scope of the in- 
fluence of various design parameters upon system performance. The linear 
approximation is proposed as a useful design guide when applied with due 
awareness of its limitations. In addition, a design criterion is presented 
by means of which it is possible to avoid designs that lead to operation 
in a "hovering" region. Such operation leads to prolonged delays in 
ball opening time. Application of the design criterion leads to a 
ball/projectile design that is optimum with respect to ball opening time. 
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INTRODUCTION 



BACKGROUND 

As early as 1858, the potential benefits of tubular projectiles were 
recognized for specific purposes. In that year, Joseph Whitworth (Whitworth 
Threads) wrote about and included an illustration of the tubular projectile 
in the section on Rifled Firearms of his Miscellaneous Papers on Mechanical 
Subjects [1]^. The projectile devised by Whitworth was hexagonal, with a 
circular hole, and he noted its particular "...effectiveness in perforating 
elastic materials which prevented them from closing up." Whitworth further 
observed that the tubular projectile, which utilized a wooden sabot, penetrated 
deeper into masonry than any with which he was acquainted. 

According to Charters and Thomas [2], the "Krnka-Hebler" projectile 
was reported in the Allgemeine Schwei tzerische Mi 1 i tarzei tung as having been 
so successful that the Ordnance Department in the United States carried out 
firing tests of caliber .30 tubular projectiles in 1894. A description of 
those test firings and an updating of the results, performed by I.E. Segal, is 
included in Ref. 2. In the 1894 tests, a vulcanized fiber sabot was used to 
seal the projectile while in the barrel. 

The purpose of the 1894 test was to compare trajectories of standard and 
tubular projectiles. This was done by comparing the vertical drop on target at 
a given range. Segal's report indicates that the drag coefficients computed 
from the 1894 results agree closely with those of Charters and Thomas [2], 
even though the latter work was reported fifty years later (1944) and utilized 
more sophisticated experimental apparatus. 



^Numbers appearing in brackets refer to the list of references. 
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RECENT DEVELOPMENTAL STUDIES 

The result obtained by independent researchers as to the relative merits 
of the tubular projectile as an alternative to the conventional round are con- 
tradictory. It seems as if there was as much research devoted to discrediting 
the tubular projectile as there was to investigating its merits. Frank and 
McLaughlin [3] have accumulated a great deal of data from various sources in an 
attempt to "objectively compare" the merits of the tubular and "conventional 
shapes". The authors concluded that the tubular has no particular advantages 
over "well designed" conventional projectile shapes. Their findings disagree 
with many of those discussed in the following sections. 

Range Tests 

Winchenbach, Daniel and Edgar [4] conducted range tests of six configurations 
of tubular projectiles and concluded that the drag coefficients were significant- 
ly lower than the standard High Explosive Incendiary (HEI) projectile of the same 
caliber. Only projectiles of the same bore size were compared. The experimental 
models were constructed from standard 20-mm ammunition by boring various hole 
sizes to obtain "...area ratios (At/A-j) of 0.7, 0.8, 0.9, and 1.0." Where A^ is 
the throat area and A^ is the inlet area of the projectile. The drag coefficients 
of the projectiles with area ratios of 0.8, 0.9, and 1.0 were less than half that 
of the projectile with area ratio of 0.7. The higher drag coefficient was attri- 
buted to the underexpanded flow at the exit resulting in high base drag. 

Range tests conducted for the concept evaluation of the 20-mm tubular pro- 
jectile for the Vulcan Weapons System [5] showed that the tubular projectile per- 
formed better than the standard HEI round, with lower drag, shorter time of flight 
for a given range (30 percent shorter at 1000 meters, 40 percent shorter at 2000 
meters) and similar dispersion characteristics on target. 

Recent tests at the Naval Weapons Center, China Lake, California of the Ball- 
Obturated Tubular Projectile (BOT) have shown a definite advantage in that the low 
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drag of the tubular projectile means slower retardation of velocity and hence 
a higher terminal kinetic energy [6]. (The BOT is the object of this investi- 
gation and will be described in detail in later sections of this report.) 

Target Impact Tests 

The higher the energy delivered to the target, the greater the damage. 

Since kinetic energy is directly proportional to the first power of mass and 
to the square of the velocity, if a tubular projectile has relatively less mass 
this maybe compensated by a higher muzzle velocity. However, a reduction in 
mass to obtain higher muzzle velocities may not be necessary, depending upon 
the imposed performance criteria. 

Target impact studies have been carried out with various projectile con- 
figurations and target types. Rhethorst, et.al . [7] conducted impact studies 
of 7.62-mm tubular projectiles on helmets. Tests showed that even with the 
same energy of impact, the tubular projectile penetrated further. 

Kitchen and Keeser [8] conducted studies for the Air Force on the impact 
effectiveness of tubular projectiles on simulated aircraft fuel cells. These 
tests were conducted with steel and depleted uranium (DU) tubular projectiles 
and standard 20-mm HEI projectiles. The projectiles were fired at double 
panels at varying degrees of obliquity. Of the forty-seven tests of the steel 
tubular projectiles which impacted the target, twenty breached the rear panel 
when fired at angles up to 70°. The DU projectiles breached the rear panel 
even at angles of 85° and fires were started in three of the five DU tests. 

The standard M56 HEI projectiles failed to breach the rear panel in each of 
five firings even though severe damage to the front panel and two fires occurred. 

Brunsvold and Kalivretenos [9] conducted a program to test the effective- 
ness of the tubular projectile against a simulated cruise missile warhead. 

Tests were conducted with 20 and 40-mm rounds. Results of those tests are 
classified, and therefore not presented here. It was noted [9] that due to the 
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improved stability characteristics of the tubular projectile, only about 
half the spin rate need be imparted by the rifling in the barrel to obtain 
the same stability as a standard projectile, with obvious implications with 
respect to wear. 

Weapons System Compatibility 

For the tubular projectile to be practical, it must be able to interface 
with existing weapons systems. Reference 5 is a report of the 1978 tests 
conducted by the Army to evaluate the tubular projectile in the Vulcan weapons 
system. 

Because the nose of the tubular projectile is flat compared to the stan- 
dard projectile (See Fig. 1), it "...did not lend itself to chambering in the 
weapon." Personnel safety required that the weapon be remotely operated and 
therefore only the surface-to-surface performance comparison was made. 

The dispersion on target improved for the tubular projectile at ranges of 
from 700 to 2000 meters. The dispersion varied on the average by only 0.2 
mil over the entire range spectrum between the tubular and conventional HEI 
projectile. However, at 2000 meters, the tubular was better by 0.4 mil. The 
tubular projectile had a 30 percent shorter flight time at 1000 meters and a 
40 percent shorter flight time at 2000 meters than the conventional HEI 
projectile. 

Results of firing of 25-mm and 30-mm tubular projectiles from the Oerlikin 
KBA and GAU-8/A respectively were reported in part in Ref. 3 Figures comparing 
Line-of-Si ght Penetration Capability vs. Range for a conventional spinner con- 
figuration based on the AR-2 shape, a finner configuration based on the FIO 
shape, and a tubular projectile were shown. The tubular projectile was shown 
to be inferior in both instances. The report does not indicate the degree of 
compatibility the tubular configuration has with the guns used. 

Internal Shock Wave Considerations 

From photographs of tubular projectiles in flight [4] and during wind 



4 



tunnel testing [7], various shock patterns have been observed at the inlet 
region and in the wake. The photographs of Figs. 2a-c show the BOT and the 
conventional projectiles in flight. The sabots used in Refs. 4 and 7 were 
of the pusher type consisting of a simple base plug. The detached bow 
shock seen in the photograph of the standard projectile (Fig. 2a) is similar 
to that of the tubular projectile shown in Fig. 2b with the tubular passage 
closed (and in Ref. 4 with the sabot attached). With no blockage of the 
internal passage (Fig. 2c) there was no detached bow shock, with only 
attached oblique shocks emanating from the lip of the projectile. This latter 
case allows the possibility of four other internal flow/shock configurations. 

(1) A normal shock standing at the entrance to the projectile; 

(2) A normal shock standing at some intermediate position in the 
channel ; 

(3) A system of oblique shocks present in the channel; 

(4) The channel is devoid of shocks and the flow is shock-free 
through the channel . 

Reference 10 contains a discussion of the shock patterns and the 
various conditions leading to their existence. 

PURPOSE OF STUDY 

The results of tests and experiments previously mentioned show, for 
the most part, definite advantages accruing to the tubular projectile over 
the conventional projectile primarily in the area of lower drag (approximately 
1/2 to 1/3 of conventional [8]) which allows a flatter trajectory and shorter 
time of flight for a given range [5]. With reference to air combat and 
anti-missile defense both from a Surface-to-Ai r and Air-to-Air standpoint, 
these features are definitely worth further study. 

The problem seen from the air platform is that of the discarding sabot 
or pusher. The solid [4, 5, 7,] or split [8] disk is effective and reliable. 
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Figure 1. 20-nm Ball Obturated Tubular Projectile (30T) 

components, and conventional 20-mm projectile 
(Photograph courtesy of NWC, China Lake) . 




Figure 2a. Standard 20-mm projectile in flight with 
detached bow shock. (Photograph courtesy 
of NWC, China Lake) . 
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Figure 2b. 30T with passage blocked resulting in 

detached bow shock. (Photograph 
courtesy of NWC, China Lake) . 




Figure 2c. EOT with no blockage in passage; only attached 
oblique shocks emanating from the lip of the 
projectile. (Photograph courtesy of NWC, China 
Lake) . 
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however, the possibility of ingestion in the aircraft engines prohibits 
their use. Rhethorst et al . , in Ref. 7 illustrates a number of sabot/ 
obturator designs for use with the tubular projectile, but these are also 
of the discarding type. They also investigated the possibility of using a 
consumable sabot [11] which would be burned up as the projectile was tran- 
siting the barrel and would be completely consumed as it exited the barrel. 

The most recent development has been the BOT which was designed at 
NWC China Lake, California [6]. The ball obturator has been bored with a 
hole the same diameter as the hole through the projectile. When loaded, 
the ball is supported inside the projectile such that the axis through the 
hole in the ball makes an angle of approximately 90® with the axis of the 
hole through the projectile (see Fig. 3). Gas pressure from the burning 
propellant holds the ball obturator fixed with respect to the projectile 
(as well as to the ball). Upon exiting the barrel, the propellant gas 
pressure is released and, in its place, forces due to aerodynamic effects 
come into play. At that time a complex force distribution acts on the ball 
to cause it to change orientation inside the projectile in such a way as to 
align the holes. Inertial forces then dominate to maintain this alignment 
so that the projectile remains fully tubular. There are no separating 
parts and the opening process is automatic. 

The purpose of the study reported herein has been to analytically predict 
the motion of the ball obturator as a function of time. In support of this 
purpose, an experimental apparatus was designed and built to simulate the 
motion of the obturator in order to provide physical evidence of the accuracy 
of the model and to gain insights into the nature of the motion. 
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Figure 3 




Sketch showing dimensions and positioning of the ball 
obturator within the spherical cavity. 
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AfWLYSIS 



DEVELOPMENT OF THE EQUATIONS OF MOTION 

The ball obturator is a rigid body of revolution with a system of 
coordinate axes fixed to the ball (body-fixed axes) having its origin at 
the mass center of thg ball (Fig. 4). These axes are designated x, y, and 
z, where z is the axis through the hole in the ball. 

To describe the motion of the ball relative to the projectile an 
inertial frame of reference is defined with axes designated X, Y and Z 
with its origin also at the mass center of the ball. For this analysis 
it is assumed that the projectile is constrained to rotate about the Z-axis. 
Because the motion of the ball with respect to the projectile is of special 
interest it should be noted that components of the ball angular velocity 
with respect to the X and Y axes are also relative to the projectile, 
whereas the component with respect to the Z-axis must account for the pro- 
jectile spin. 

The motion may be described by Euler's Modified Equations of Motion 
for a rigid body of revolution about a fixed point. These are given in 
Thomson [12] as: 



Aw 


+ (C-A)w w = M 


(la) 


1 


2 3 1 


Aw 


+ (A-C)w w = M 


(lb) 


2 


1 3 2 




Cw = M 

3 3 


(Ic) 



where 

A = Mass moments of inertia about principal axes 

perpendicular to the z-axis. 

C = Mass moment of inertia about the principal 

axis through the hole (z-axis). 
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y 



Figure 4. niustration of coordinate systems and symbols used in the 
analysi s . 



V> Z 
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0 ) , ci) , 0 ) = Angular accelerations about the body-fixed axes 

1 2 3 

X, y, z respectively. 

CO , oj , oj = Angular velocities about the body-fixed axes. 

1 2 3 

M , M , M = External moments acting about the body-fixed 

1 2 3 

axes . 



The mass moments of inertia for the particular case at hand 



(sphere of radius R with concentric hole of radius 


r) are given by: 


A [%(l-r2)5/2 (3^12 + 2) + (1 - r2)5/2 


(2) 


C = (3^2 + 2) 


(3) 


where r = r/R. 




The symbol X is introduced here to denote the 


relative difference 


between the two moments of inertia. That is: 




, C-A 




^ - -5- 




From Eqs. (2) and (3) it is seen that 




5p 


(4) 


f^+4 


and since 0 - r - 1, X 1 ies in the range 0 - X - 


1. 


With X so defined, Eqs. (1) may be rewritten: 




CO + Xco CO = M /A 


(la) 


1 2 3 1 


CO - Xco CO = M /A 


(lb) 


2 13 2 


u = M /C 


(Ic) 



The orientation of the ball at any time may be described by 
three angles; 0, (p and if). These are known as Euler angles and their 
relationships to the coordinate axes already described may be seen in 
Fig. 4. The components of angular velocity co , co and co in terms of 

^ 2 3 

the Euler angles are given as [12]: 
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oa 

1 


= ’4)sin9sinf{) + 9cosi 


(5a) 


0) 

2 


= ii'Sin9coS(i) - 9sinc{) 


(5b) 


0) 

0 


= 4> + ii»cos9 


(5c) 



By differentiating Eqs. (5) with respect to time, the angular 
accelerations, w , w , and u may be obtained as 

1 2 3 

0 )^ = ii)sin9sin4) + ii»(9cos9sin(j> + |sin9cos4>) + 9coS(j) (6a) 

-9*sin(}) 

w = ii)Sin9cos4> + ii)(9cos9cos4> - |sin9sin4>) - 9sin^) (6b) 

2 

-9|cos(t) 

w = (p + i{>cos9 - '|)9sin9 (6c) 

3 

Combination of Eqs. (1), (5), and (6) yields the Euler equations 
in terms of the Euler angles: 

ii»sin9sin(j) + iii(9cos9sin^> + isin9coS(j)) + 9coS(J> - (7a) 

9|sin({) + X(ii»sin9cos<j) - 9sin4>) (ii»cos9 + |) = M^A 

if)Sin9coS(}) + i(9cos9cos<^ - |sin9sin<j>) - 9sin<J) - (7b) 

9|cos<}) - X(il)sin9sin(i) + 9coS(l)) (iiicos9 + |) = M /A 

2 

ti>cos9 - :J)9sin9 + = M^/C (7c) 

If Eq. (7a) is multiplied by sin({) and Eq. (7b) by cos(}), and the 
results are then added, we obtain^ 

iljsin9 = (M /A)sin(}) + (M /A)coS({) + 9 [(X + 1)| + (X-l)iiicos9] 

1 2 



■^The overbar notation is introduced here to denote dimensional time- 
dependent quantities. These will subsequently be non-dimensional ized 
through the use of the initial projectile spin rate. 
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Similarly, if the multiplication is reversed and the results are 
subtracted, we have 

9 = (M /A)cost{) - (M /A)sincf) - ipsinG [(A + 1) + AiJicosG] 

1 2 

Simply rewriting Eq. (7c) we have 

n “ ~ n 

(j) = R /C + i|i9sin0 - i|;cos9 

3 

These are the basic governing equations for the motion of the ball 
in terms of the Euler angles. They are second order ordinary 
differential equations, nonlinear, and strongly coupled. Before 
proceeding further, we shall introduce the convenience of non- 
dimensional ization of the time-dependent variable. This is ac- 
complished through the use of the initial rate of spin of the 
projectile. The equations of motion are unchanged in form except 
for the terms involving the external moments. 

Thus, with the definitions 

M = R /AaJ„„^ M = M /Aw„„^ M = R /Co) ^ 

1 1 PO 2 2 PO 3 3 PO 

we have 

4;sin9 = Misincf) + Macosc}) + 9 [(A + 1) | + (A-1 )iIjcos9] (8) 

9 = MiCOS4> - MaSincf) - iiisin9 [(A + 1)(J + A4icos9] (9) 

> = Ms + )i)9sin9 - i|)COs9 (10) 

— •• n 

where ip = 'P/cOpoi ^ > etc. 

The integration of Eqs. (8) - (10) proceeds from a condition 
in which the ball is fixed to the projectile and spinning with it. 
The initial values of the precession angle and the spin angle <p 



14 



are arbitrary and are set to zero for convenience. The appropriate 
initial conditions are, therefore: 



0(0) = e 



0 



0 ( 0 ) = 0 



( 11 ) 



4 >( 0 ) = 0 



Ho) = 0 



i{^(0) = 0 



t^(O) = 1 



The central problem yet to be addressed is that of the specification 
of the external moments M , fi , and M or, in non-dimensional form, 

: 2 3 

M , M , and M . 



EXPRESSION OF THE EXTERNAL MOMENTS 

The interaction of the ball v/ith the cavity in the spinning 
projectile is dominated by torques arising from the relative motion 
between the two bodies. These torques will be opposite to the direction 
of the relative velocity so that we may begin by specifying that: 



where 1^ = = unit vector in the direction of the relative 



We have considered torques due to both fluid shear and sliding friction. 
Torque Due to Fluid Shear, M^ 

Nakabayashi [13] has recently reported the results of a series of 
experiments to determine the fluid shear acting in the gap between tv/o 
concentric rotating spheres. Although the details of the flow are 



2 



M = -1 M 

~ r 



( 12 ) 



angular velocity 
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extremely complex [14], the results of Nakabayashi show that the 
data are well correlated by: 

C = 8ir(l + e)/(3eRe) for laminar flow (13) 

m 

and = 0.053[1 + (7/4)e]/{e[l - (3/2)e]Re° ' for turbulent flow, 
where = M^/(pR®co^^) , e = (h/R)«l, and the Reynolds number is given 
as : 



Re = R^ai^/v. 

From [13] it is also possible to estimate a transition Reynolds 
number: 

Re^ = 70e'^'® 

In our experiments (both real and numerical) we have found that relative 
angular velocities in the range of interest are well below those 
necessary to stimulate turbulent flow in air so that Eq. (13) is 
adequate for our purposes. When Eq. (13) is cast in a suitable form we 
have: 



= 8ir(l + e)viR^o)y(3e) (14) 

Torque Due to Sliding Friction, 

In the case of sliding friction the magnitude of the applied torque 
may be written 



Me = h r^ 
s s n m 



(15) 



where F^ is the normal component of the resultant of all forces acting 
upon the ball and r^ is the moment arm extending from the point of con- 
tact and directed normal to the vector of relative angular velocity. 

Since the point of contact must be on the surface of the ball, 
the position vector of this point, P, must have a magnitude of the 
ball radius R. Figure 5 illustrates the geometry for cases to be dis- 
cussed in which multiple points of contact are distributed along a line. 

If Y and 3 are the location angles of a point of contact, as illustrated 
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X 




Y 



Figure 5. Sketch illustrating the location of points of contact between 
ball and projectile. 
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in Fig. 5, then the position vector may be written 
P = R(e singcosy + e sin3siny + e cos3) 

-V a/ 



where e , e , and e are unit vectors for the inertial coordinates. 






The vector of relative angular velocity may be written, 



0 ) = eco +eo) +eco 

, x rx ry 7z rz 

and, since P • co = Rco cosn> we have 

a, -v' 

n = cos"^ (o) sinBcosy + o) sinBsiny + co cos3)] (16) 

Oj^ ia ry I L 



The angle n is that subtended by the position vector P and 

'\j 

the relative angular velocity vector co so that the moment arm r^^ 

'\j 



is given by 

r^ = Rsinn (17) 

and from Eq. (16) we obtain; 




With the point of contact specified the moment arm is therefore 
determined and, together with the normal force F^, leads to the torque 
due to sliding friction, Eq. (15). 

At present it is impossible to exactly specify the vector of the 
forces applied to the ball. We have therefore examined two special 
cases that are relevant to the investigation. These are: 

1. Motion under a gravity load. 

2. Motion under an axial aerodynamic load. 

Gravity load . In this case the load is the weight, W, of the 
ball and is applied directly downward. The point of contact is at the 
bottom of the cavity so that 3 = tt/o and y = tt (that is, P = -Re ). 

The moment arm is given by 

r„,n = - (w 

nig rx r 
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( 19 ) 



and the magnitude of the moment is 

Which is applied opposite to the relative velocity vector in accordance 
with Eq. (12). 

Axial aerodynamic load . The aerodynamic load acting upon the 
ball will be a complex function of the ball orientation. During the 
opening instants, however, when the ball blocks the tubular projectile, 
the resultant of the pressure forces on the ball will act along the axis 
of the projectile. This interval is of importance, of course, since it 
will constitute a significant part of the ball-opening delay. At lower 
values of 0, when there is significant flow through the ball, the specif- 
ication of the net aerodynamic force will require further detailed analysis. 
It is important to note further that when the ball is partially open there 
may be a significant moment acting upon the ball, as well as a centrally 
directed force. An understanding of the nature of the aerodynamic force 
under partial ly-open conditions is a major goal of continuing research. 

For the case of an axially oriented applied force, the ball will 
be forced against the aft rim of the hole in the projectile. The exact 
number and location of the contact points is difficult to predict and 
will, in fact, depend upon projectile manufacturing variations. However, 
any effect of arbitrariness in contact point specification is "blurred" 
by the high-speed rotary motion and errors resulting should be random 
and therefore of little influence on the overall motions. This hypo- 
thesis has been verified by computer experiments. In selecting the number 
of contact points and their location, we have assumed the first stable 
position, that is, three contact points. In addition, we have assumed 
that these points are symmetrically dispersed along the line of contact. 
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The axial aerodynamic force, F, , is evenly balanced among the 

a 

contact points so that at each point the normal force is given by 

F = ^ cosB 

" 3 

and, since the contact points lie along the rim of the projectile 
hole, 3 = sin~^ r/R. (see Fig. 5). The moment due to sliding friction 
in this case is given by 

( 20 ) 

where 3 = sin’^(r/R) and the three values of r^^ are given by Eq. (18) 
with the angle y set at 0, 2tt/3, and 47t/3. 

The aerodynamic force most appropriate for use in Eq. (20) will 
be that due to a pressure distribution over the ball and behind a normal 
shock standing at the projectile entrance, (see sketch below) 




nrTT 



CO 



The net force on the ball will be 

F = •iTr^(P - P ) 
a 2 1 

Since the actual base pressure P^ is much less than P 2 and, in any case, 
is unknown, we assume ambient pressure here. That is, P^ = P^. The face 
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Pa will lie somewhere between the static and the stagnation pressure 
existing behind the shock. The difference is small (about 7%) in the 
range of Mach numbers of interest ( = 3 to 5) and will likely 

approach the static pressure because of spillage behind the shock and 
recirculation within the cavity ahead of the ball; that is, there will 
be little or no pressure recovery in the cavity. The assumption of 
static pressure is adopted for these reasons as well as the fact that 
this is a conservative assumption in terms of ball opening times (higher 
pressures lead to shorter opening times). Thus, the approximate ex- 
pression for the pressure P 2 is 



Equations (20) and (21) provide the necessary information for determining 
the torque due to an axial aerodynamic load and this torque is, again, 
applied in a direction opposite to that of the relative velocity vector. 
Consideration of the Relative Angular Velocity \fector 

The three body-fixed components of the relative velocity vector are 
given by see Eqs. (5) : 

0 ) ^ = (ii» - Wp)sin6sin(f) + 0cos(f) 

CO = (tIi - CO )sin0coS(i) - OsimJ) (22) 

2 r p 

“3^ ^ (’i' ■ “p)cOS0 + 

According to Eq. (12), then, the components of the external moments 
may be written: 

H = - -T - co„)s in0sin(i) + 0cos*l 

1 %o “r ^ 

M^ = ' ^ 2 ^ [(lii - cop)sin0cos(}) - 0sin(})] (23) 




and for the aerodynamic force. 




( 21 ) 



po 



M 
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The terms involving these moments in the equations of motion may 



now be written 



-M 



[(’!' - a)p)sin9] 



M sin(|) + M co5i> = 



1 2 




M coscj) - M sintf) = 



(24) 



1 2 



M = as in Eq. (23) . 



The value of R in these expressions are given according to the table 
below: 



It is appropriate to combine the fluid shear moment with either 
of the sliding friction moments. In most cases analyzed, however, 
fluid shear moments are several orders of magnitude smaller than 
those due to sliding friction. Since the gravity and axial aerodynamic 
loads are each dependent upon specified points of contact, it is not 
appropriate to combine them without determining a point of contact that 
is correct for the combination. In cases of practical interest, however, 
the aerodynamic load is either absent or, if it is present, it completely 
dominates any effects due to gravity. Consider, for instance, a sea level 
launch under the following conditions: 
r = 4.7 mm 
R = 7.9 mm 



flo = 3.0 

Under the approximation that the moment arms are about equal to the sphere 



Source of moment 
Fluid Shear, R^ 
Sliding Friction 
Gravity only, Rsg 



Applicable Equations 
(14) 



Axial Aerodynamic, R_, 

S a 



(19) 

( 20 , 21 ) 
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radius for both gravity and aerodynamic loads, the ratio of moments 
for the two cases is 






For the stated conditions with a steel ball, W =.08 N (about 0.4 lb.) 
and F = 65 N (about 15 lb.) so that the aerodynamic effect is clearly 

d 

dominant. 

It should be noted, in closing this section, that the main 
weakness in the analysis is due to uncertainties in specifying the 
aerodynamic load. If this load (force vector or centrally applied 
force plus moment) is known, the methods given to predict its affect 
upon the motion are quite general. 
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EXACT SOLUTIOMS FOR TYPICAL SYSTEM PARAMETERS 



Equations (8) - (10), subject to the initial conditions 
[Eq. (11)] have been solved for both the gravity-only case and the 
case of an axial aerodynamic load. The gravity-only case is of special 
interest because it provides theoretical predictions for the experimen- 
tal program (described in a subsequent section). The aerodynamic load 
is, of course, the practical case for the tubular projectile. 

The appropriate equations were coded for solution on the digital 
computer. The general flow of the solution method is as follows: 

1. Set initial conditions and specify system parameters. 

2. Calculate necessary constants (ball weight, moments of inertia, 
flight level pressure, etc.). 

3. Calculate oj , oo , co from current Euler angles [Eqs. (5)]. 

1 2 3 

4. Calculate body-fixed components of relative angular velocity 

[Eqs. (22)]. • 

5. Convert relative angular velocity components to inertial axes 

(o) , CO , CO ) . 

V rx’ ry’ rz' 

6. Determine moment arms and external moments as required (see 
table, p. 22). 

7. Determine moment terms, for use in equation of motion [Eqs. (24)]. 

8. Determine \p, 9, and <j> [Eqs. (8) - (10)]. 

9. Integrate for new values of 4>, 9, |, and ip, 9, and 4>. 

10. Return to step 3, until 9 is sufficiently close to zero. 

The method of integration employed (step 9 above) was a fourth-order 
Runge-Kutta scheme with a step-size automatically adjusted to control 
numerical precision. 



24 



Consideration of Projectile Angular Velocity . 

Two system constraints have been investigated in the analysis. 

In the first case we have invoked the constraint of constant pro- 
jectile spin rate and, in the second case, we have required a con- 
stant system (ball plus projectile) kinetic energy. Since the in- 
ertia of the ball is small relative to that of the projectile (a 
factor of about 20) the decrease in projectile spin rate necessary 
to maintain constant system kinetic energy during the ball motion 
is rather small (typically about 2%) and the two cases are not sig- 
nificantly different. The computer solutions shown here (so-called 
"exact" solutions) are for the case of constant system kinetic energy. 

Gravity Load . 

The system parameters for this case were set to correspond to 
the conditions of the experiments to be described later. These para- 
meters are summarized in Table I. 

Figures 6 and 7 illustrate the nature of the motion predicted by 
the exact solutions for a representative set of the physical parameters. 

In Fig. 6 the angle of tilt (9) relative to the cavity spin axis is 
shown as a function of time for various rates of spin. At relatively 
high rates of spin the motion appears to be quite linear with a nearly 
constant rate of nutation (0) developed early in the motion. At lower 
rates of spin, however, there is an initial period during which the mo- 
tion is nearly stable, followed by a rapid decrease in the tilt angle. 

The time required to reach a given value of 9 is therefore a function of 
the angle itself and is illustrative of the complex interactions between 
the inertial tendencies of the ball and the dissipative torques due to 
friction. Figure 7 is included to show the variations in the Euler angle 
rates for a typical case. When the ball is released its rate of precession 
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TABLE I 



Ball Characteristics and Experimental Conditions - Gravity Load 



Ball Designation 


Mod 0 


Mod I 





Physical 

Characteristics 



R, m 


7.899 X 10"^ 


7.925 X 10'^ 


r, m 


4.722 X 10"^ 


1 .984 X 10’^ 


C, Nms^ 


3.126 X 10"^ 


4.029 X 10'^ 


X 


0.410 


0.077 


P, kg/m^ 


7.67 X 10^ 


7.75 X 10^ 



Experimental 

Conditions 

1.917 X 10’^ 
0.30 

6.43 X 10"^ 
89.5 
14.50 

2.18 X 10'^ 



^Empirically determined, see text. 
2 

Computed from Eq. ( 39 ). 



£ 

9 

0 

9* 



Ns/m^ 
(note 1 ) 



deg 
deg 

M, Nm (note 2) 



1.917 X 10 
0.51 

6.43 X 10 

87.25 

36.71 



-5 



-5 



2.08 X 10' 
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06*0 09*0 




Figure 6. Exact solutions for 0 vs t. Gravity load only; ball Mod 0. 
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.00 




Figure 7, Euler angle rates (non-dimensional) for typical 

case. Gravity load only, ball Mod 0, = 400 s’^ 
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(ij)) increases above its initial non-dimensional value of unity and, 
concurrently , a retrograde spin is developed. For the conditions of 
Fig. 7 (wpQ= 400 s"M» the mean value of nutation rate is fairly constant 
with acceleration rates that are initially large and rapidly damped. 
Additional results for case of gravity load are discussed in con- 
juntion with the linear approximation and experimental results. 

Axial Aerodynamic Load . 

Figures 8 and 9 present similar results for the motion under an 
aerodynamic load. The system parameters for these figures are those 
of the "baseline case" listed in Table II. 

The nature of the response is seen to be similar to that of the 
gravity-only case. In particular, the existence of a hovering mode is 
apparent at low projectile spin rates. This phenomenon points up some 
design limits and is the subject of discussion in a subsequent section. 
Sensitivity studies were conducted for the aerodynamic load and these 
also are discussed later. 
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TABLE II 



Ranges of Design Variables - Aerodynamic Load 



Design Variable 


Minimum 


Basel ine 


Maximum 


Figure No.^ 





R, mm 
A, Nms^ 


basel ine 


7.90 

2.24x10"’ 


39 5 

7.60x10"‘‘ 


24 


r/R 

A, Nms^ 
X 


0.400 7 

3.20x10" 

0.19 


0.598 

2.24x10"’ 

0.41 


0.800 
1. 00x10"’ 
0.69 


25 


P , kg/m^ 

A, Nms^ 


3.10x10^ 

8.96x10"® 


7.75x10® 

2.24x10"’ 


1.24x10'* 

3.58x10"’ 


26 




0.71 


0.855 


1.000 


27 


M 

00 


3.00 


3.50 


5.00 


28 


atm 


0.1 


1.0 


basel i ne 


29 


po 


6.91x10® 


1.20x10“^ 


l.OOxlO"* 


8 


0 , deg 


70 


85 


basel ine 


8, 10 



^ Figures cited are those in which the given design variable is 
investigated over a range of values. 
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t» ms 
(b) 



Figure 8. Exact solutions for 9 vs t; baseline projectile, (a) High 
projectile spin rates; (b) low projectile spin rates. 
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Figure 9. Euler angle rates_for a typical case. Baseline projectile, 
cOpQ = 1.20 xlO** s ^ (115,000 rpm). 
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LINEAR APPROXIMATI ONS 



From Figs. 8a and 8b it will be noted that there is a marked change 
in the nature of the ball motion at a projectile spin rate of about 
7300 s'^ (70,000 rpm). At rates above this value (Fig. 8a) the ball motion 
appears to be quite linear; that is, the response of the ball to the friction- 
al torques approximates one of a second-order linear system responding to a 
combined ramp and step input (both of which are negative). Such a response 
is given by 



+ exp(-^t) {Rsin[(l-^2)*5^ ^ _ osin[(l-c^)^t + <J>p)]} /(l-;;^)^ 

where <5^ = = tan’^[(l-?^)*^/^] 

The parameters of the response, to be derived below, are given by; 

C = M/2, D = -(A/2)sin26^, and R = MD/(A + 1) 

The method leading to the linear approximation is a perturbation upon the 
equilibrium solution to Eqs. (8) - (10) under torque-free conditions. 

DEVELOPMENT OF THE PARAMETERS OF THE LINEAR RESPONSE 

Equation (25) satisfies the initial conditions on 6 and is the solution 
to 



From the equations of motion for 6 and the associated relation for external 
moment [Eqs. (9) and (24)], a comparison of the coefficient of the tilt ve- 
locity, 8 , leads to the damping coefficient: 



0 = + D + R(t - 2c)... 



(25) 



e + 2c0 + 6 = 6^ + D + Rt 
0 



(26) 



C = M/2 



(27) 



where M = M/Ato o) . 

po r 

With this definition of c , combination of Eqs. (9) and (26) leads 
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to 

(9 - 9^) - (D + Rt) = ii)sin9[(A + I)! + Xiiicos9] (28) 

The response specified in Eq. (25) indicates that the function described 
in Eq. (28) is oscillatory in nature with the mean value decaying exponenti- 
ally from an initial value of -D to an asymptotic final value of -2R^. That 
is, the mean value (about which the oscillations occur) is given by; 

[(9 - 9j - (D + = -2Rz - (D - 2R?)exp(-i;t) 

In the limiting case of vanishingly small external moments this expression 
simplifies to a straight-line response with slope R so that an evaluation 
of this rate of decay should emerge from an investigation of the solution of 
the equations of motion for the case of zero external moments. 

If the external moments are removed from Eqs. (8) - (10), a solution 
of the following form is obtained: 

(X + 1)4)^ + Acos9^ = 0, 9^ = 0, ’I’q = 1 (29) 

The subscript (^) is used here to indicate quantities pertaining to this 
solution. Perturbations to this solution, designated by subscript {^), are 
defined such that 

M'i^p 

<j> = 4 >q + M4>p (30) 

9=9 + M9„ 

0 p 

The base quantities (subscript o) must satisfy Eq. (29) and the perturbed 
quantities (no subscript) must satisfy the initial conditions, Eq. (11). 

The appropriate initial conditions for the perturbations are, therefore 
'I'p(O) = 0 

lp(0) = 4 q/M = -Xcos9^/M(A + 1) (31) 

6p(0) = 9p(0) = 0 
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When Eqs. (29) and (30) are inserted into Eqs. (8) and (10), the fol- 
lowing approximate relationships are obtained: 

Mi(» sin6 = T, - M6 cos9 + O(M^) 

^p 0 1 p 0 ' ' 

= T» + M0 sin6 - Mii) cos6 + O(M^) 

p 2 p 0 p 0 

where the torque terms and are obtained from Eqs. (24), (29), and (30) 



T.J = -M(MiiipSine^) = 0(N|2) 

= -M^^/(X + 1) + O(M^) 

Neglecting terms of order and smaller, we have 

ijj sin6 = -0 COS0 
p 0 p 0 

$ = -(^ /(X + 1 ) + 0 CSC0 

p 0 P 0 

The first integrals of this set yield 



4 ' 

o 

4 ) 



P 

P 



= -0 cot9 
P 0 

= -(})^t/(X+ 1) + 0pCSC0^ 




(32) 



or, from Eqs. (30) 

= 1 - M0„cot0„ 

p 0 

(33) 

^ = -M(j)^t/(X + 1) + M0pCsc0Q 

When Eqs. (29) and (30) are substituted into Eq. (28) we have, to the first 
order : 

Msin0^[(X + l)4>p + Xii/pCOS0^ - X0pSin0^] = M0p - (D + Rt) 

and with Eqs. (32) this becomes 

(X + 1 )sin0^[M|^t/(X + 1) + (|i^] = D + Rt 



It follows that: 

D = (X + 1)$ sin0 = -(X/2)sin20 (34) 

'00 0 

R = MD/(X + 1) 
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In order to evaluate the non-dimensional moment parameter, M, it is 
necessary to develop approximate relationships for and for the moment 
arm r^ to be used in the previously derived expressions. Under the con- 
straints of the linearization described here (M«l and lD|«l), the 
relative angular velocity is approximated by [Eqs. (22)]: 

“r " ^3r ~ *^o " “ ^cos9^/(X + 1) 

The principal component of the relative angular velocity is therefore 
along the spin axis of the ball and, since the ball is initially near 
9^ = ir/2 (required for |D|«1) the components along the inertial axes are, 
approximately, 

“rx^“r “ ='*"*0 
“rz/“r ” “ 

Since is a function of time (even though 4/^ = constant = 1 in the 

linearization), suitable mean values for the moment arm relationships must 

be obtained. In the case of a gravity-only load we have 

- (oo /to )^]^= R|cosip I 
mg rx r ' o' 

the absolute value being necessary to preserve the positive nature of the 
moment arm. For the mean value, 

1 

■•mg “ j ''mg 

we have 

- 1 2r 

~ / R|cosi|j I di|/ = — / cosi|j dip 

mg 2ir /q ' ^o ' ^o tt 

and 

r^g = (2 /tt)R (36) 

For the case of aerodynamic loads, we have from Eq. (18) 
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''ma ^ " sin^3sin^(iij^ +y)]^ 

where the selected values of locate the contact points and B = sin'^(r/R). 

In the mean, then 

7 R ^ ^ ^ 1 

^ma = — /q - sin^Bsin^(ij;^ +y)r d.].^ 
and 

^ma " 

where E(B) is the Complete Elliptic Integral of the second kind. 

With these values for the moment arms it remains to evaluate the line- 
arized forms of the external moments. In the case of fluid friction we 
have, from Eq . (14) , 

= [8tt( 1 + £ )yR^to^/(i; )]/Au^(i)pQ 

or 

= 8tt( 1 +£)yRV(^Ao3p^) (38) 

Note that no linearization is required for this non-dimensional moment para- 
meter. For the case of gravity-only loads, from Eq. (19) 

= [y2WR(2/Tr)]/Ad)^a3pQ 

or 

Msg = 2 usWR(A + D/TrAAWpQ^cose^ (39) 

For the aerodynamic load Eq. (37) gives, for three contact points symmetrically 
di sposed , 

3 

= 3(2/tt)R E(b) 

and from Eq. (20) 

Msa = y5FgR(2/7T)E(B)cosB/A(I)^a3pQ 

or 

M.. = 2y F R(A + 1 )E( B)cosB/ttAAw ^cos0^ (40) 

Sa S a pO 0 
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The fluid friction term, may be added to either the gravity term, 
or the aerodynamic term, M , , to obtain M and the related parameters 

Sa 

[^, R, and D from Eqs. (27) and (34)]. 

RESULTS COMPUTED FROM THE LINEAR APPROXIMATIONS 

Figure 10 illustrates the linear approximation in comparison with the 
exact solution for two initial values of 0 under fluid friction and gravity 
load. (As mentioned previously, the effects of fluid friction are generally 
negligible. This observation also applies in the linearization.) Figure 11 
gives similar results of the case of an axial aerodynamic load on the ball. 

It is seen from both figures that good approximations to the exact solution 
are obtained from the linearization. The approximation departs from the 
exact solution as time passes and the ball angle decreases. This is to be 
expected since the terminal phase of the ball motion is not adequately 
modeled in the linearization. Nevertheless, the linear model provides valu- 
able insights into the nature of the motion and its most influential parameters. 
It should be noted again that at the lower angles, when an observer would 
see a sizable opening through the projectile when sighting along its axis, 
both theoretical models are subject to the possibility of considerable error. 

The is especially true in the aerodynamic case because of uncertainties in 
the load description at partial ly-open conditions. 

The results of the linear approximation are discussed further in con- 
nection with the experimental program to be described in the next section. 
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Figure 10. Comparison of approximate solutions with exact solutions. 
Gravity load only, ball Mod 0, = 400s~^. 
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Figure 11. Comparison of approximate solutions with exact solutions. 
Baseline projectile, = 1.20x10'* s"^. 
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EXPERIMENTATION 



EXPERIMENTAL APPARATUS 

An experimental system was designed to simulate the spinning BOT 
under the influence of a gravity load. The apparatus consisted of a 
compressed-ai r-dri ven spin-up rig, optical-timing mechanism, air manifold 
and associated piping, tubing and electronics. 

The obturator was fabricated from a standard 5/8-in. diameter chrome- 
steel (52100) bearing ball. The ball was annealed to allow machining and 
then bored along a diametrical axis. The bored ball was then mounted in a 
bakelite metalographic specimen mount. The mounted ball was placed in a 
milling machine and a flat was machined in the bakelite at a specified angle 
relative to the axis through the hole in the ball. This flat was then used 
as the polishing plane for metalographic specimen preparation. A small flat 
spot was polished on the obturator to provide a highly reflective surface at 
a known orientation relative to the z-axis of the obturator (Fig. 12), 

After removal from the bakelite mount, the obturator was placed in a 
three-piece, lucite housing (Fig, 13). The mating ends of the two hollow in- 
ner cylinders were each machined with a 5/8-inch spherical end mill to a depth 
of approximately 5/16-inch. When mated, a spherical cavity was formed to 
accomodate the obturator. The third cylinder was press fit over the others 
after the obturator was inserted to insure alignment of the inner pieces and 
provide rigidity (Fig, 14). The ends of this lucite assembly were then press 
fit into aluminum end-pieces similar to those in Fig. 14. These end-pieces 
served as the shaft for the bearings and one also served as the prime mover 
(bucket wheel) for the apparatus. The shaft rotated in two ball bearings 
mounted in aluminum pillow blocks aligned on a rigid pedestal. The prime mover 
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Figure 12. Polished spot orientation on obturator and relation- 
ship to 9 . 

0 



42 




Figure 13. Obturator and components of lucite obturator 
housing with mating ends of inner cylinders 
machined to house the obturator. 




Figure 14 . Assembled obturator housing with aluminum 
end pieces . 
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was a bucket wheel machined from a solid aluminum disk (Fig. 15). 

The bucket wheel was driven by compressed air supplied from an install- 
ed system through an air filter to a manifold and then through one of two 
Model-10 Kendall pressure regulators. From the regulator, the air passed 
through a flexible tube to a tee, each leg of which supplied a nozzle. These 
nozzles were mounted opposite one another on the pillow block in such a way 
as to allow the air jets to impinge upon the bucket wheel to cause rotation 
(Fig. 16). The speed of rotation was sensed by a Bentley Nevada Proximitor, 

Model 3100N, which was mounted above the bucket wheel to detect the passage 
of each point on the wheel. The pulses generated by the sensor were counted, 
averaged and displayed as a frequency by a Monsanto Programmable Counter-Timer 
Model HOB. 

While the spin-up rig was being brought up to the desired speed, an 
air jet(subsequently referred to as the holding jet) from a nozzle mounted 
rigidly on the pedestal at the end opposite the bucket wheel held the obturator 
fixed to the spinning rig. The holding jet passed through a hole in the 
plunger of a Rocker Solenoid, R.S. No. 10-207. This hole was aligned with 
the hole through the "projectile" to allow the jet to impinge upon the obturator 
(Fig. 17). When the appropriate switch was activated, standard 115VAC was applie 
to the Rocker solenoid and removed from a normally-closed ASCO Solenoid Valve 
in the air supply line to the hoi ding- jet. The plunger retracted from the 
position shown in Fig. 17 and the solenoid valve closed. The retraction of the 
plunger performed three functions. It first caused a pulse to be generated 
by a proximitor mounted next to the nozzle. This pulse started the timer 
function of a second Model HOB Counter-Timer. The plunger, by misaligning its 
hole with that of the nozzle, blocked the holding jet. This removed the re- 
straining influence of the holding jet on the obturator and prevented any air 
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loading" position. 
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Figure 16. Spin-up rig showing bucket-wheel with air je 
orientation and location of Proximitor above 
bucket wheel. 







47 



Figure 17. Holding-jet nozzle (1), Proximitor (2) and Solenoid 
plunger (3) extended prior to obturator release, 
(hole through plunger is aligned with the end of 



remaining in the supply line from impinging on the obturator. And lastly, 
in the fully retracted position, the machined and polished end of the plung- 
er was positioned opposite the axis of the apparatus (Fig. 18). When the 
obturator, now free to move relative to the "projectile", nutated through 
enough of an angle, the beam of a SPECTRA PHYSICS MODEL 132 LASER, MODEL NO. 
3187, passed through the projectile and was reflected by the polished end of 
the plunger onto a light sensitive diode (Fig. 19). The pulse emitted by the 
diode triggered the stop channel on the timer and the elapsed time was dis- 
played. Figure 20 shows the entire experimental set-up with air flasks in the 
background. 

EXPERIMENTAL PROCEDURE 

The obturator was positioned in the spin-up rig in such a way that the 
laser beam was reflected off the polished spot back to a target mounted on the 
laser (Figs. 20 and 21). The center of the target is the location of the beam 
and the obturator was adjusted to place the reflected spot on the center of 
the target. Thus the orientation of the z-axis through the obturator was 
known relative to the z-axis (laser beam) through the projectile. This is the 

initial value 0 . 

0 

Once the alignment had been checked and the obturator positioned and held 
at the known 0^, flow to the turbine nozzles was controlled in such a way that 
the desired spin rate was achieved. The turbine flow regulator was then adjust- 
ed so as to maintain the spin rate at ±5 Hz on the digital display or approxi- 
mately ±1.75 rad/sec. Once the desired rate was achieved and noted, and the 
timer checked and reset if necessary, the solenoid switch (Fig. 17) was acti- 
vated. The obturator was released and began to move relative to the spin-up 
rig. When the angle 0 reached the value at which the laser beam could pass 
through the obturator hole [0=0* = sin'^r/R)], the beam struck the polished 
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Figure 19. Light sensitive diode housing ( 1 ) illuminated by 
reflected laser beam from the polished end of th 
solenoid plunger. 
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Figure 21. Laser reflection on the target from the polished 
spot on the obturator. 
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plunger and was reflected into the light sensitive diode housing. The 
reflected beam striking the diode caused the diode to emit a pulse which 
turned off the timer gate. The elapsed time, from plunger retraction to 
release to the obturator nutating to 9*, was displayed on the timer dis- 
play to the 0.0001 second. This time, t*, was recorded. 

EXPERIMENTAL DETERMINATION OF THE COEFFICIENT OF SLIDING FRICTION, u . 

s 

From Eq. (25) it will be noted that for large values of Cico t the time 

po 

elapsed for nutation from 9^ to 9* is given by: 

t* = -[(9^ - 9*) + D]/(Ro3p^) + 2?/03p^ 

Without the influence of fluid shear (which is negligible in this case), 

Eqs. (34), (35) and (39) lead to: 

t* -»■ 7rAcOpQ[(9^ - 9*) - (X/2 )sin29Q]/(2p^WRsin9^) 

as 2;/(Op^- 0 

At relatively high spin rates, therefore, the elapsed time t* is expected 
to vary directly with cOp^. Such being the case, the coefficient of sliding 
friction, may be estimated from the slope of data for t* vs 9* at high 
val ues of CO 

po 

DISCUSSION OF THE EXPERIMENTAL RESULTS 

Figure 22 shows the results of experiments in which the angular velocity 
was varied within the range of capability of the experimental apparatus. 

Two balls were tested, each with different hole sizes (see Table I) but 
otherwise essentially identical. Also shown in Fig. 22 are the theoretical 
predictions given by the exact solution and by the asymptotic form of the 
linear approximation, Eq. (41). This latter expression has been matched to 
the data at high spin rates in order to obtain empirical estimates of the 
coefficient of sliding friction, p^. 
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Figure 22. Presentation of experimental results showing comparisons 
with exact and approximate solutions. Gravity load only. 
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The experiments verify the theoretical prediction of the high-speed 
asymptote and the existence of a minimum in the elapsed time for ball -open- 
ing, with respect to projectile spin rate. The discrepancies between theory 
and experiment at lower spin rates (where the linear approximation is not ex- 
pected to be valid) are thought to be mainly due to inadequacies in the ex- 
perimental apparatus. At these low spin rates there is a prolonged period 
of "hover" of the ball near its initial position, followed by a rapid de- 
crease in tilt angle (see Fig. 8). During these initial instants of time 
any external disturbances will tend to cause a premature conclusion of the 
hovering period and it is this effect that is thought to have led to the low 
measured values of t* at the lower values of co„„. The effect is more pro- 
nounced in the Mod I ball which, due to its small value of X, is considerably 
less prone to align itself, particularly at the lower values of cavity spin 
rate. 

The experiments were also limited in precision due to the interaction 
of the holding jet with the solenoid-actuated plunger. This difficulty 
contributed to the scatter of data and limited the maximum usuable value of 
because of the high jet pressures required to hold the ball in position 
at high spi n rates . 

The empirically determined value of the coefficient of sliding friction, 

0.30 for the Mod I ball, appears to be reasonable in comparison with the value 

of 0.35 reported for laminated plastic on steel [15]. The larger value of = 

0.51 for the Mod 0 ball is thought to be due to the influence of the relatively 

large hole which, in the high-speed rotating motion, will have an effect 

similar to that of roughness. Within the precision of the experiments there 

was no detectable influence of speed upon and the assumption of a constant 

value for this parameter appears to have been justified, at least for the 

higher values of o)„^. 

= po 
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SENSITIVITY TESTS AND DESIGN GUIDES FOR 
AERODYNAMIC LOADING 



Figure 11 shows a comparison of the linear approximation with the 
exact solution for the motion of the baseline model (Table 2) at the 
nominal spin rate of 115,000 RPM. Two initial ball tilt angles are also 
shown and it may be seen that the linear approximation yields results that 
are quite acceptable for initial angles near 90°. The values of M, R, and 
D are given on the figure for reference purposes. The usefulness and limita- 
tions of the approximate method are also illustrated in Fig. 23 where the 
initial opening time is given as a function of projectile spin rate. As 
in the experimental program this time, designated as t*, is the time required 
for a straight-through path to open along the projectile centerline. The 
ball tilt angle at this instant is given by 

9* = sin'M>f'/R) 

and the corresponding time is a useful reference quantity for comparison of 
various obturator designs. (In addition, the aerodynamic force model used 
in this analysis becomes highly suspect at values of 0 below 9* since the 
simple shock structure assumed here will be invalid if there is any signifi- 
cant flow through the projectile). 

Figure 23 indicates that the linear approximation approaches the exact 
solution at large values of This is to be expected since M decreases 

quadratical ly with cOp^. The asymptotic behavior of t* at high Up^ is a 
useful design insight and is given by 
= Aco„«[( 9„ - 9*) + D]/Msin0^ 

°° PO 0 0 

and for the baseline case this reduces to t* - 4.0 x 10"‘*oj ms 

” V po 

(cOpQ in s The sharp minimum in t* illustrated in Fig. 23 suggests some 
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Comparison of approximate solutions with exact solutions. 
Baseline projectile, = 1. 20x10"* s‘^ 
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caution in the design of these devices and also provides encouragement re- 
garding the existence of criteria leading to an optimum design. The low 
spin-rate behavior shown in Fig. 23 is outside the linear range of operation 
and is due to the delaying mechanism of sliding friction at relatively large 
values of M (large friction, low inertial imbalance due to a small hole in 
the ball, and/or low projectile spin rate). 

Under these conditions the ball "hovers" at its initial value and con- 
siderable time elapses before a ball orientation is achieved in which friction 
is in a direction to aid the alignment of the hole with the projectile axis. 
This behavior has been observed in previous solutions (see Fig. 8), and in 
the experimental program (see Fig. 22). Design guidelines necessary to avoid 
the hovering mode are discussed below, following a discussion of the results 
of the sensitivity analysis. 

SENSITIVITY ANALYSIS 

The design variables affecting the obturator performance are those con- 
tained in the normalized torque determined by Eqs. (18) - (21). The influence 
of the time-varying relative angular velocity and moment arms in these ex- 
pressions may be qualitatively evaluated by referring to the linear approxi- 
mation so that the relevant design variables are [Eq. 40] X, y^, F^, g, 0^, 

R, and cOp^. Since X and g are geometric parameters depending upon the radii 
of the ball and its hole; and the minor moment of inertia. A, depends upon 
these as well as the density of the ball material, this list may also be 
written: r, R, p^ y^, M^, P^, and Here the aerodynamic force has been 

replaced by means of Eq. (21) and the projectile hole is assumed to be of the 
same size as that of the ball. In considering the coefficient y^, it must be 
remembered that it depends upon the two materials (ball and projectile) and 
is therefore dependent upon p. We have studied independent variation in y^, 
however, since its value is subject to considerable uncertainty. Experiments 
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have shown that the value of should be approximately 50% larger than 
that published for dry friction between smooth surfaces. Table 2 indicates 
the baseline value for these design variables as well as the maxima and minima 
that determines the range of sensitivity analysis. In each case (except 
C0po» which has been previously discussed) the values of the variables were 
chosen such that equal increments were provided and the baseline condition 
was included in the range. Single parameter variations are shown here al- 
though work is continuing using the techniques of numerical optimization. 

The computer data are presented as 9 vs t and the interested reader 
may construct cross-plots such as that of Fig. 23 to evaluate the effects 
of these design variables upon the time required to reach a given angle, such 
as 8*. Figures 24-29 illustrate the relative importance of the selected 
design variables. It is important to note that in these figures the range of 
system performance has been restricted to the quasi -linear range previously 
discussed where the ball response is typically of the type shown in Fig. 8 
and above the minimum in Fig. 23. Outside of this range (in the "hovering" 
state) the trends illustrated may be reversed. Attempts to reduce ball re- 
sponse time by the design improvements suggested by Figs. 24-29 may push the 
performance beyond the minimum illustrated in Fig. 23, with resulting "negative 
improvements". In many cases the implications of Figs. 24-29 are self-evident 
but to assist in physical interpretation a few explanatory remarks are offered 
below. 

Geometry effects . 

Both the ball radius, R, and the radius of the ball hole, r, were varied 
independently. Figure 24 shows the effects of increasing the ball radius while 
maintaining a constant value of r/R. The variation in this case is essentially 
one of scale-up which covers a range of projectile diameters from 20mm to 100mm. 
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Figure 24. Effect of ball size on ball response. 
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It may be seen that larger systems respond more slowly. Increases in applied 
torques are not sufficient to offset the effects of the fifth-power relation- 
ship between ball inertia and radius. Figure 24 indicates that opening of the 
ball obturator in large projectiles may be unacceptably slow unless compensat- 
ing design actions are taken. As will be seen in subsequent figures, such 
actions might include increases in ball hole radius and reduction in the den- 
sity of the ball material. 

Increasing the size of the hole in the ball, with all other parameters 
held constant, reduces the inertia of the ball and increases the imbalance 
betv/een the moments of inertia about the two principal axes. Both effects 
are conducive to rapid ball response and this is illustrated in Fig. 25. 

In addition, the initial opening angle, 0*, is directly increased by an 
increase in r/R with a corresponding decrease in the elapsed time. 

Material effects . 

Though the bal 1 -obturator must be of sufficient strength to plug the 
projectile during launch, it does not contribute to the structural integri- 
ty of the projectile in flight. The ball material may therefore be open 
to some design variation and this possibility is investigated in Fig. 26 
where it is seen that the time required to reach a given ball angle is 
approximately proportional to the density of the ball material, p A 
change in ball material from Carbon steel (p 'v<7,750 kg/m^) to an Aluminum 
alloy (p 'V'3,100 kg/m^) would lead to a halving of the response time. This 
effect also leads to an opportunity of linearly tailoring the ball response 
time by means of its density--a feature that may be particularly useful in 
applications of the system as a hydromechanical switch. 

Coefficient of sliding friction . 

. As has been mentioned, the correct value of is subject to consider- 
able uncertainty. Figure 7 shows the impact of an uncertainty range of 

about ±20% in which, in terms of t*, is about ±15%. In applications for 
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Figure 25. Effect of ball hole size on ball response. 
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Figure 26. Effect of density of ball material on ball response. 
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Figure. 27, Effect of coefficient of sliding friction on ball response. 
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which margins of error such as this are unacceptable, tests of prototype 
configurations may be necessary to determine The range tested in 
Fig. 27 represents a best estimate for dry friction of steel on steel 
for the range of r/R covered in Fig. 25. Systems using dissimilar materi- 
als will require appropriate adjustments in u^. In addition, it may be 
seen that some roughening of the interfacial surfaces may not be detri- 
mental to system performance while, on the other hand, intential polish- 
ing of these surfaces may lead to slower ball action. 

Launch conditions . 

Projectile launch Mach number and altitude have been separately in- 
vestigated with the results shown in Figs. 28 and 29. The increase in 
wave drag associated with higher Mach numbers are actually beneficial in 
terms of ball response, because of the corresponding increases in the 
external moments that aid the motion of the ball within the projectile. 

On the other hand, increases in launch altitude lead to increases in 
ball response times. In proceeding from sea-level to a flight-level 
pressure of 0.325 atm (about 30,000 ft in the standard atmosphere), for 
example, a threefold increase in t* is predicted for the baseline design. 

DESIGN GUIDES 

The linear approximation previously developed represents a signifi- 
cant simplification to the computational tasks required for the exact 
solutions. As such, this simplified approach is recommended for design 
and "first-cut" calculations. As has been seen, however, there is a 
minimum in elapsed time (see Fig. 23, for example) which represents a 
demarkation of the limit of validity of the linear approximation as well 
as a point of optimum design for the ball obturator. A means for esti- 
mating this optimal point is developed below. 
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Figure 28. Effect of launch Mach number on ball response. 
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Figure 29 . Effect of launch level pressure on ball response. 
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Under the initial conditions of the problem at hand ('i»=l , 4>=0, 9=0, 

9=9^) the equations of motion [Eqs. (8) - (10)] may be written: 

ip = 9^cot9„ 

0 0 

S = -9 cos29 /sin9 
^000 

9 = -M9^ - (A/2)sin29^ 

Here the term M9^ is retained since it may be significant for large values 

of M even though 9^ is initially zero. In addition, in the initial instants 

of motion, the relative motion of the two bodies is largely due to nutation. 

That is, (D == -9 , and 
r 0 

M9 = -M/Ao) ^ 
po 

Thus an approximation to the initial acceleration 9 is given by 

9 = (M/Aw + D = constant 
po 

In cases in which hovering occurs this expression will be approximately 
zero and an indication of a hovering design is therefore given by 



M 

Aw ^ 
po 



-0 



and the design guide to avoid hovering is given by 
M 



Aw "nDl 
po ' ' 



< 1 



In the initial instants, the applied moment due to aerodynamic loading 
may be conservatively estimated as = y^F^RcosB so that the design 

Scl S a 



criterion is given as 



H = 



2 y F^RcosB 

s q < 



AAw ^sin29 
po 0 



< 1 



(42) 



where, as usual, B = sin~^(r/R). Note that in this expression as A 0 



or 9^-»- tt/ 2 the factor H grows beyond bound. The first situation occurs 
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if there is no hole in the ball (A = C) and the second case represents the 
stable initial condition. Equation (42) indicates a hovering condition in 
both cases, in accordance with expectations. 

As an example, consider the baseline case (Table II). For this 
design and Mach number, - 93 N (about 21 lb). Inserting the other 
constants from Table II we have 
H = 0.735 

and the prediction is that hovering is avoided. If the same design is con- 
sidered but Eq. (42) is solved for in the critical case (H = 1) we 
have C 0 p^(cri tical ) = 10,000 s"^ or about 98,000 rpm. Comparison of these 
calculations with the results indicated in Figs. 8 and 23 shows a very 
good agreement. In particular. Fig. 23 shows that if the equality is 
used in Eq. (42) a very good estimate of the optimum design is obtained. 
Thus, for a given worst-case launch condition (maximum F, due to maximum 

a 

Mach number and minimum altitude) and initial projectile spin rate (low- 
est value for worst case) the optimum ball design (X, A, B, 9^) may 
be deduced from Eq. (42). It would be wise, of course, to apply some 
safety factor to this optimum design and, in any case, a violation of 
the inequality in Eq. (42) would lead to hovering - a failure in the 
operation of the bal 1-obturator . To repeat, if the inequality is satis- 
fied in Eq. (42) the linear approximation [Eqs. (25) and (40)] may be 
safely appl ied. 
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CONCLUSIONS 



The gyrodynamic motion of a ball obturator within a spinning tubular 
projectile has been investigated from both theoretical and experimental 
points of view. The theoretical developments have included both exact and 
approximate formulations and comparison of these two theories has led to 
the definition of the range of validity of the approximate method. 

The experimental program was conducted under conditions in which the 
main source of external moment was sliding friction due to the weight of 
the ball. The results of these tests have substantiated the validity of 
both the exact and approximate theories and have led to several insights re- 
garding the motion. In addition, a semi-empirical method of measuring the 
coefficient of sliding friction evolved from the experimental program. 

Two distinct regions of motion have been identified. The first of 
these - the hovering mode - is characterized by a metastable initial ori- 
entation for prolonged periods of time. This undesirable behavior may 
be avoided if design guides developed within this study are followed. 

The second region of motion occurs when the hovering region is avoid- 
ed and in this mode, the response of the ball is approximately linear. The 
analysis has shown that minimum ball -opening times may be expected for ball 
designs and projectile launch conditions that are near (but not beyond) the 
onset of hovering motion. 

The analytical model has been used to conduct sensitivity tests to de- 
termine the relative influence of design parameters for ball designs that are 
in the proper (non-hovering) range. In general, ball response is more rapid 
for larger hole sizes, lighter material, greater sliding friction, higher 
launch Mach number, and lower launch altitudes. Small obturating balls will 
respond more quickly than large ones, all else being equal. The theoretical 
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predictions depend directly upon the nature and magnitude of the aero- 
dynamic load acting upon the ball in flight. Uncertainties regarding the 
theoretical description of this load are the main source of doubt re- 
garding the precision of the predictive models presented herein. 
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LIST OF SYMBOLS 



A Mass moment of inertia about minor axis 

C Mass moment of inertia about major axis (axis of hole in ball) 

D Step input associated with linear model 

e Unit vector 

/Si 

F Force 

h Gap width separating ball and cavity (fluid friction model) 

H Design criterion, see Eq. (42) 

k Ratio of specific heats 

M External moment (also used to denote Mach number) 

P Pressure 

r Radius of hole in ball 

R Radius of ball 

R Slope of ramp input associated with linear model 

r Moment arm 

m 

t Time 

Greek Letters 

6 Angle used in locating contact points ( = 6* = sin~V/R ) 

Y Angle used in locating contact points (see Fig. 5) 

e Dimensionless gap width ( = h/R ) 
t; Damping coefficient associated with linear model 

n Angle used in locating contact points (see Fig. 5) 

A (C - A)/A 

y Molecular viscosity of fluid in gap 
p Density of ball material 
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LIST OF SYMBOLS 






<P,'I>,Q 

to 

“p 

1,2,3 

a 

f 

g 

n 

0 

r 

s 



x,y ,2 



00 



* 



Greek Letters (continued) 

Phase angle associated with linear approximation (see Eq. 25) 

Euler angles 

Angular velocity 

Spin rate of projectile 

Subscripts 

Refer to body fixed axes 
Aerodynamic quantity 
Fluid quantity 
Refers to gravity load 
Normal component 

Initial quantity (also refers to solutions of torque-free motion) 

Denotes relative angular velocity 
Associated with sliding friction 
Refer to inertial (projectile-fixed) axes 
Vector quantity 
Conditions in free stream 

Superscripts 

Dimensional quantity - absence infers non-dimensional ized time-dependent 
quantity. E.g. 0 = ^(^pq) • 

Refers to condition of first opening of 1 ine-of-sight along projectile 
axis and through hole in ball obturator 
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